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Microwave response of an NS ring coupled to a superconducting resonator 
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A long phase coherent normal (N) wire between superconductors (S) is characterized by a dense 
phase dependent Andreev spectrum . We probe this spectrum in a high frequency phase biased 
configuration, by coupling an NS ring to a multimode superconducting resonator. We detect a 
dc flux and frequency dependent response whose dissipative and non dissipative components are 
related by a simple Debye relaxation law with a characteristic time of the order of the diffusion 
time through the N part of the ring. The flux dependence exhibits h/2e periodic oscillations with a 
large harmonics content at temperatures where the Josephson current is purely sinusoidal. This is 
explained considering that the populations of the Andreev levels are frozen on the time-scale of the 
experiments. 



PACS numbers: 

Most properties of a non superconducting (normal) 
metal connected to two superconductors (an SNS junc- 
tion) can be seen as resulting from the Andreev states 
(AS) in the normal metal, and the occupation of those 
states. Andreev states are correlated electron-hole eigen- 
states in the normal metal, determined by boundary con- 
ditions imposed by the superconducting banks. In a long 
diffusive normal metal (of length L grater than the su- 
perconducting coherence length S,s), the AS spectrum 
is a quasi continuum of levels with a small energy gap 
Eg [1, 2]. This so-called minigap depends solely on L 
and the diffusion constant D, via the Thouless energy 
Eth ~ 'hlT-o where tu = L^/D is the diffusion time 
along the junction. Eg is fully modulated by the phase 
difference (p between the superconducting order param- 
eters on both sides. Eg(ip) is maximal at </? = with 
Eg{0) ~ S.lETh and goes linearly to zero at ip — tt, 
as was recently measured by scanning tunneling spec- 
troscopy [3]. The phase dependent Josephson current 
I J (if) at equilibrium sums the contributions of each An- 
dreev state of energy e„, via i„ = -^g", the current 
carried by the level n of occupation factor p„. 



ri'P) =^Pn{<^n{'P))in{f) 
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As was measured recently by Hall magnetometry [6], 
the phase dependence of the supercurrent is non sinu- 
soidal at low temperature, but turns sinusoidal at T > 
Eg, with an amplitude that decreases roughly exponen- 
tially with temperature on the scale of Eg . 

Whereas the equilibrium properties of the proximity 
effect (PE) in long SNS junction are well understood 
theoretically and experimentally, its dynamics is a more 
complex and yet unresolved issue. In particular the cur- 
rent response of the system to a time dependent phase 
which affects both energy levels and their populations 
still needs to be determined. We expect at least the exis- 
tence of two time constants, one related to the relaxation 



of the populations to a change of the energy levels, an- 
other one related to the dynamics of the Andreev levels, 
and possibly also a contribution due to the generation of 
quasiparticles. Which of the electron-electron, electron- 
phonon, diffusion or dephasing times are relevant in these 
processes? Experimentally, there are many ways in which 
to impose an out of equilibrium situation. Because of the 
Josephson relation linking the voltage to the time deriva- 
tive of the phase, a voltage bias across the junction causes 
a time dependent phase, and the I{V) curve is a probe of 
the dynamics of the PE. Irradiating the junction with an 
rf field, which produces Shapiro steps in the I{V) curve, 
is another. Both the critical current in current biased 
experiments [4, 5] and the complete Josephson current 
phase relation [6] have been shown to be very sensitive 
to RF excitation at frequencies of the order of or larger 
than I/tjj. In addition, the electron-phonon scattering 
rate has also been shown to set the frequency scale for 
the hysteresis in current biased experiments [5]. 

In contrast to all these experiments performed in 
strongly non linear regimes, in this Letter we present 
a linear response experiment in which we measure both 
the non dissipative (inphase), x' and dissipative (out of 
phase), x" current response of an SN ring. We measure 
these with imposing an oscillating phase difference, at 
several RF frequencies, as a function of the dc supercon- 
ducting phase difference. Our results identify a current 
relaxation mechanism with a relaxation time of the or- 
der of the diffusion time. The temperature is greater than 
the minigap, so that the dc current phase relation is sinu- 
soidal. Nevertheless, we find that the dissipative and the 
non dissipative components contain several harmonics. 
This shows on the one hand that at high frequency, the 
non dissipative response is not simply the flux derivative 
of the supercurrent, and also that dissipation is enhanced 
when the minigap closes. 

The experiment consists in inductively coupling an NS 




FIG. 1: Electron micrograph of the NS ring inserted in the 
niobium resonator. 



ring to a multimode superconducting resonator operat- 
ing between 300 MHZ and 6 GHZ. In the linear response 
regime the dc flux dependences of x' and x" are deduced 
from the flux induced variations of the resonator's res- 
onances frequency and quality factor Q. A similar tech- 
nique based on single mode resonators was already used 
for the investigation of short Joscphson junctions inbed- 
ded in superconducting rings [8, 9]. 

The resonator consists of two parallel superconducting 
Nb meander lines (2/L(m wide, L^ = 20 cm long, and 4^m 
apart, see Fig. 1) patterned on a sapphire substrate [12]. 
One of these lines is weakly coupled to a RF generator 
via a small on chip capacitance whose value is adjusted 
in order to preserve the high Q of the resonator(80 000 at 
20 mK), the other one is grounded. The resonance condi- 
tions are Lji — ^, where the mode index n is an integer, 
and A the electromagnetic wavelength. The fundamental 
resonance frequency (n — 1) is of the order of 360 MHz. 
The resonator is enclosed in a gold plated stainless-steel 
box, shielding it from electromagnetic noise, and cooled 
down to mK temperature. The high Q factors allows us 
to detect extremely small variations: -4- ~ SQ^^ ^ 10^^. 
It therefore provides very accurate ac impedance mea- 
surements of mesoscopic objects. This technique was 
previously employed in contactless measurements of the 
response of Aharonov Bohm rings [7] . A single NS ring 
gives more signal than 10^ normal rings, for two reasons: 
First, the supercurrent in an NS ring is g times larger 
than the persistent current in a normal ring of the same 
size as the N wire (g is the conductance of the normal 
wire in 2e^//i units). Second, the length of the S part can 
be adjusted to optimize the inductive coupling between 
the ring and the resonator. 

The N part of the NS ring is a mesoscopic Au wire 
(99.9999% purity, 1.2 x 0.38 fim'^, 50 nm thick, normal 
state resistance {Rn = ^/Gn = 0.52 fi) corresponding 
to 5 = 2 X 10"*), positionned between the two meander 
lines by e-beam lithography, and thermally evaporated 
onto the saphire substrate. We estimate the Thouless 
energy for this sample to be Eth = 90 ± 10 mK. The 
S part of the ring is constructed by conecting the Au 
wire to one of the resonator lines via W nanowires de- 
posited using a focused Ga ion beam (FIB) to decom- 
pose a tungstene carbonyle vapor. The superconducting 
transition temperature of the W wires is 4K, their crit- 




FIG. 2: Flux dependence of the real x' and imaginary suscep- 
tibility x" for different temperatures and resonance frequen- 
cies (left panel, f=365MHz, right panel T=0.67K). The am- 
plitude of the Josephson currents calculated are 7j(0.67 K) = 
18 nA and Ij{lK) = 6 ^A . The data are corrected from ge- 
ometrical inductance effects, which induce corrections which 
are sizable at 0.67 and 0.55 K but negligible at IK . Circles: 
fit with Xd{T)/{l -\- iujt) see expressions 3 and 5 calculated 
from Usadel equation [17] at temperatures corresponding to 
ksT = 5, 6 and QErh giving the best agreement with ex- 
perimental data. The amplitude of the theoretical curves has 
been rescaled by a factor 0.7. Note that this experiment only 
measures the flux dependence of x, all the curves have been 
arbitrarily shifted along the vertical axis. 



ical current is 1 mA, and critical field is greater than 7 
T [10]. This technique provides a good interface qual- 
ity, thanks to a slight etching step prior to deposition 
of W. We have checked that SNS junctions fabricated 
using this technique exhibit supercurrents and Shapiro 
steps comparable to long SNS junctions made with more 
conventional techniques [11]. The sample is thus an ac- 
squid embedded in the resonator. The dc supercon- 
ducting phase difference (p at the boundaries of the N 
wire is imposed by a magnetic flux $dc created by a 
magnetic field applied perpendicularly to the ring plane: 
(f — ~2TT^iic/^o where 4>o = h/2e is the superconduct- 
ing flux quantum. In addition, an ac flux S^i^cos{ujt) 
is generated by the ac current in the resonator. At low 
enough frequency the current in the SNS ring should fol- 
low adiabatically the oscillating flux, and the ac response 
should be entirely in phase, given by the flux derivative 
of the Josephson current 9/j($dc)/9$ , also called the 
inverse kinetic inductance. But at high frequency the 
ac current SI^j in the ring should have both an in-phase 
and an out-of-phase response to the ac flux : 6Iuj{^dc) ~ 
X'(w, <^dc)S'^u,cos{ujt)+x"{^, ^dc)S^u;Sin{ujt). Here x = 
x' -I- ix" is the complex susceptibility of the ring related 
to its complex impedance Z by x — ii-^/Z. The aim of 
the experiment is to determine xi^dc) of the ring at the 
successive resonances of the resonator by measuring the 



flux dependences of the resonance frequencies and quality 
factor which are related to xifn) via: 
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C = 0.15 fiH is the total inductance of the resonator, 
M ^ 5 pH the mutual inductance between the ring 



and the two lines of the resonator and Lg 



15 pH 



the ring's geometrical inductance. The coefficient fc„ 
accounts for the spatial dependence of the ac field am- 
plitude at frequency /„. The limit Lg — yields the 
well known formulas where Sf/f and 5{\/Q) are respec- 
tively proportional to x' ^-i^d x"- But since the geo- 
metrical inductance of the ring is finite, the internal dc 
flux $i„t differs from the applied flux '^ext and reads 
'^int — '^ext + Lglj{^int)- This Screening effect leads to 
Eq.2 and to an hysteretic behavior at low temperature 
when the parameter /3 = 2TrLgIc{T)/^o is larger than 
1. (The critical current Ic is the maximum amplitude of 
/j($)). Because of this, it is not possible to access the 
whole current/phase relation below 600 mK [12]. 

In the following we present the dependences of x upon 
the internal flux through the ring. Our main findings are 
summarized in Fig. 2 showing the flux oscillations of x' 
and x" for different temperatures and resonator modes. 
Each curve corresponds to an average of 30 to 100 mag- 
netic field scans. Their periodicity corresponds, as ex- 
pected, to $0 through the ring. The first surprise is the 
observation, beside the non dissipative response x'i^)^ 
of a larger dissipative response x"($), over the entire 
frequency range investigated (365 MHz to 3 GHz). In 
addition, x'(<&) and x"{^) are not harmonic functions of 
flux even though the dc Josephson current is sinusoidal 
in this temperature range. It is clear in particular that 
x'(<I') is not equal to xji^) = 9/j($)/9$. It is also inter- 
esting to compare the flux variation amplitudes 6x' , Sx" 
to the inverse kinetic inductance Xj{T) ~ 27r/c(r)/$o 
and to Xjv = wGat . Whereas we find that 5x! is smaller 
than Xj(^)7 "^x" is niuch larger than x'n- We now turn 
to the frequency dependence x{'^) ™ order to extract the 
relevant relaxation times. The amplitudes 5x' and 5x" 
are plotted in Fig. 3 for the successive resonator eigen fre- 
quencies /„ = w„/27r, at two temperatures. We find that 
the frequency dependence follows a simple Debye relax- 
ation law, (5x(w) = (5x(0)/(l + jcjt), with two adjustable 
parameters: 5x{Q), of the order of xj, and the relaxation 
time T equal to 0.6 ± 0.2 ns , which is 7 times the diffu- 
sion time. We could not detect significant temperature 
dependence of this time between 0.5 and IK. 

So far we have discussed exclusively the linear response 
regime (microwave excitation powers between 10~^^ and 
lO"'^^ W). Figure 4 presents non linear effects in a larger 
NS ring with the same N wire. They show up at greater 



powers correponding to an amplitude of the ac flux of the 
order of $0 through the ring: (5(1/(5) ($) features sharp 
peaks at odd multiples of $o/2, for which the minigap 
vanishes. The peaks widen with increasing microwave 
power. We attribute this extra dissipation to microwave 
induced Andreev pair breaking, which should occur pref- 
erentially in the range of dc flux where the minigap is 
the smallest. Raising temperature or frequency decreases 
the excitation amplitude threshold for the appearance of 
these peaks. The in phase response is also modified, with 
a striking increase of 5/($) in the whole flux range and 
a sign change of the flux dependence in the vicinity of 
$0/2. Non linearities also show up at lower excitation 
powers when temperature is increased. This behavior 
recalls the strong modification of the dc current/phase 
relation under microwave irradiation [6, 19]. 

We now discuss possible explanations of the linear- 
response data. One cause of dissipation in SNS junc- 
tions is the relaxation of the population of Andreev lev- 
els which explains non equilibrium effects in voltage bi- 
ased configurations such as fractional Shapiro steps [13]. 
This relaxation is characterized by the inelastic time Tin 
which is the shortest between the electron-electron and 
the electron-phonon scattering times. It yields a fre- 
quency dependent contribution to the response function 
X proportional to the sum over the whole spectrum of 
the square of the single level current [14]. 
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where F{ip, T) = -E„ 
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can be written in the con- 



tinuous spectrum limit in terms of the spectral current 
J(e) and the density of states p{e) of the Andreev levels 
as: F{ip,T) = ^ J'^{ip,e)/ [kBTp{e))]de. The high fre- 
quency limit ofxin is: 



Xd{f,T) = lim Xtn{(^) 



dip 



-F{^,T) (4) 



which can be also written Xdi'PjT) — E„ P«a 
corresponding to time independent "frozen" popula- 
tions of Andreev states. We have calculated F{ip, T) 
[15, 17] using Usadel equations [16]. In the limit 
where fc^T ^ Exh it is approximativcly given by 
F{ip, T) = {2TrGNETh/ekBT){-Tr + mod (tt + p, 27r) - 
sgn(sin((p)) sin {ip/2)/tt) sin(i^). It is is nearly tt periodic 
with a sharp anomaly at ip = tt related to the closing 
of the minigap. This leads to an important harmonics 
content of xif) at finite frequency just as found exper- 
imentally (see Fig. 2) and discussed below. However the 
experimental flux dependence of Sx" is similar to Sx'iip) 
and can not be reproduced within this model. More- 
over the relaxation of the Andreev levels population is 
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FIG. 3: Frequency dependence of S'x and S-^' for two differ- 
ent temperatures. Continuous lines: fits according to Debye 
relaxation laws with r = 0.6 ns. Inset: frequency dependence 
of the ratio 5x!' /S'x! for several temperatures. The linear de- 
pendence in 27r/ confirms the validity of the Debye relaxation 
fits. Circles: T=1K, crosses: T=0.67K 
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FIG. 4: Non linear regime. Flux dependent dissipative and 
non-dissipative responses for different excitation microwave 
powers at T = Q.TK. This data was taken for the same N 
Au wire inbedded in a larger ring than the previous measure- 
ments. 



and r expected to be of the order of td- We show in 
Fig. 2 that the experimental flux dependence of 5'x' and 
5x" at 0.55, 0.67 and IK and 2 different frequencies are 
well reproduced by Eq.(5) with r = G.6ns as determined 
previously and Xd{v) calculated from expression (4). We 
have used a unique rescaling factor of the order of 0.7 
for all the data which can be attributed to errors in the 
estimation of the mutual coupling between the ring and 
the resonator. Recent results [22] on the finite frequency 
response of a SINIS junction where the normal meso- 
scopic wire is isolated from the electrodes by an insulat- 
ing barrier lead to a frequency dependent response with 
similar flux dependences for x' and x", as observed here 
and a characteristic time given by td- The extension of 
this work to our experimental configuration with highly 
transmitting N/S interfaces is in progress [17]. Finally 
pair breaking induced by residual magnetic impurities 
could also be important and needs to be considered. 

In conclusion we have found that the high frequency 
current response at GHz frequencies of NS rings is dra- 
matically different from the flux derivative of the Joseph- 
son current. The harmonics content of the dc flux de- 
pendence is well understood as a result of the freezing 
of the populations of the Andreev states which cannot 
follow the time dependent flux. The physical origin of 
the nanosecond time scale responsible for the observed 
dissipative response and the associated high frequency 
decrease of the in phase response seems to be related to 
the diffusion time through the normal part of the ring. 

We acknowledge help of A.Kasumov and F.Fortuna for 
fabrication using the FIB as well as useful discussions 
with R. Deblock, M. Aprili, B. Reulet and L. Glazman. 



expected to be governed by the inelastic electron scatter- 
ing time [20] identical to the Nyquist dephasing time [21] 
in a finite size quasi ID wire. ( Low-energy electrons in 
the N wire are totally Andreev-reflected from the super- 
conducting interfaces). In the limit corresponding to our 
experimental situation where T ^ gEr/kB — lOOOiiT, 
r,„(r) = T^{T) = gn/ksT which is of the order of 100 
ns i.e. much longer that the time scales (ns range) in- 
vestigated in this experiment. The characteristic time r 
entering in the fits in Fig. 3 must therefore describe a 
different and faster dynamics than the equilibration of 
the population of Andreev levels described above. The 
breakdown of the adiabatic approximation when Andreev 
levels do not follow the phase oscillations is expected to 
occur when w is larger than the inverse Ginzburg Lan- 
dau relaxation time estimated to be tgl = [t^'^ /^)td in 
a wire of finite length [4, 18]. This justifies to fit the 
experimental results at a; 3> 1/Ti„ with. 



X(w, (/?, T) = Xdi'P, T)/{1 + iujT). 
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